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Coherently enhanced radiation reaction effects in 
laser- vacuum acceleration of electron bunches 

Abstract 

The effects of coherently enhanced radiation reaction on the motion of subwavelength electron 
bunches in interaction with intense laser pulses are analyzed. The radiation reaction force behaves as 
a radiation pressure in the laser beam direction, combined with a viscous force in the perpendicular 
direction. Due to Coulomb expansion of the electron bunch, coherent radiation reaction effects only 
occur in the initial stage of the laser-bunch interaction while the bunch is still smaller than the 
wavelength. It is shown that this initial stage can have observable effects on the trajectory of the 
bunch. By scaling the system to larger bunch charges, these effects may be increased to such an extent 
that they can suppress the radial instability normally found in ponderomotive acceleration schemes, 
thereby enabling the full potential of laser-vacuum electron bunch acceleration to GeV energies. 



PACS numbers: 

41.75.Jv (Laser-driven acceleration); 41.75.Lx (Other advanced accelerator concepts); 41.75.Ht (Rela- 
tivistic electron and positron beams); 52.59.Sa (Space-charge-dominated beams); 52.50.Jm (Plasma pro- 
duction and heating by laser beams; 45.20.df (Momentum conservation). 
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1 Introduction 



The recent availability of ultra-intense laser pulses has led to increased efforts to develop novel com- 
pact acceleration sch e mes f o r electron beams. These include wakefield a ccelerators (jTaiima fc Dawson , 



f979; 



Leemans et al 



I Maksimchuk et al 



2001 



Wang et al 



2006 



2000: 



Malka et al 



Borghesi et al 



2008 



Schlenvoigt et al 



2008), thin foil irradiation schemes 



2006) and laser- vacuum acceleration concepts (jStupakov fc Zolotorev 
2008). Most studies so far have concentrated on the dynamics of the individual elec- 
trons in the considered electron bunches under the action of the externally applied electromagnetic fields, 
while consideration of the self-interaction of the bunches is limited to space-charge effects. Self-interaction 
due to the radiation emitted by the bunch itself is often ignored. Neglect of this collective radiation re- 
action is well-justified if the electron density is not too high and the electron bunch size exceeds the laser 
wavelength used. 

In the last years, however, the interaction of intense laser light with nanometer to micro meter sized atomic 



clust e rs of near solid state dens ity has become the subject of thorough investigation (jSaalmann et al. 



2006 




Krainov & Smirnov 




2002 




ted from such clusters ( 


Shao et al. 



2002). This has led to the observat i on of large numbers of electrons emit 



1996; 


Springate et al. 


2003; 


Fennel et al. 



2007 



Chen et 



2002) 



In particular, under suitable conditions the production of dense, attosecond ele ctron bunches from 



laser-irradiated clus t ers has been observe d both numerically (jLisevkina et al 



(Fennel et al 



2007 



Fukuda et al 



2010f) and experimentally 



2007). In view of these developments, collective radiation reaction 



effects may become an important factor in the dynami cs of dense ele ctron bunches and should be taken 



19571 ) that these effects can even domi- 



into account. It was realized already many years ago (jVeksler , 
nate the bunch dynamics and may be exploited as a collective acceleration mechanism. More recently, 
the ultra-intense field regime in which radiation reaction effects become dominant has been analyzed 



(Bulanov et al 



sidered 



Il'in et al 



2004) and the accelerat i on of plasma sheets assisted by radiation reaction has been con- 



2001 



Kulagin et al 



20041 ) in the context of current technological possibilities. 
In the present paper, we consider the effect of collective radiation reaction on a dense sub-wavelength 
electron bunch in interaction with a laser pulse. It is shown that the interaction can be modeled by a 



particularly simple picture of radiation pressure in the direction of the laser beam on a coherently en- 



hanced effective cross section, in combination with a viscous force perpendicular to the laser propagation 
direction. Thus collective radiation reaction may be exploited in acceleration schemes as an additional 
accelerating force as well as a stabilizing force in the transverse direction. 



Classically, the net radiative effect of the interaction of charged particles with an electromagnetic wave 
is described by Thomson scattering. This is the production of secondary dipole radiation by charges 
that oscillate due to the electric field of an incident radiation wave, which can be seen as scattering of 
part of the incident radiation by the charges. In case of an electron bunch of charge q = Ne and radius 
R — i?o much smaller than the wavelength A of the incident radiation, the N electrons in the bunch 
will essentially perform an identical oscillation, yielding coherently amplified secondary radiation that is 
identical to the radiation a point charge of magnitude q would produce. Hence the dynamics of the bunch 
as a whole can be modeled by such a point charge, as long as R <C A is fulfilled. The classical Thomson 
cross section, defined as the ratio of scattered power to incident intensity I in the non-relativ i stic li mit, is 



2007). The 



<jt = N 2 (8ir/3)rQ, with tq = e 2 /(47reoro e c 2 ) the classical electron radius (jLandau &: Lifshitz . 
factor TV 2 shows the coherent enhancement of the cross section compared to that of a single electron. Since 
directed momentum of the incident radiation is scattered into dipole radiation with zero net momentum, 
momentum is transferred to the electron bunch at a rate <ttI/c = Ft- This means that the effective force 
per electron Ft/N is proportional to N, and can therefore become significant for dense bunches. It is 
instructive to compare Ft to the ponderomotive force Fp = —Ne 2 VI/(2eom e <xo 2 ) in a laser pulse, which 
is conventionally used to accelerate charges in laser- vacuum acceleration schemes. Using a typical laser 
pulse with A = 1 /im and 100 fs pulse length, the ratio Ft/Fp becomes of the order of unity for a bunch 
of N — 10 6 electrons. Assuming a laser intensity of I = 1 ■ 10 19 W/cm 2 , this corresponds to a force of 
Ft = 22 mN, which is equivalent to an electric field of Fp/(Ne) = 0.14 TV/m in the laser propagation 
direction. This field already compares to the accelerating fields produced in wakefield accelerators, and 
increases to well above 1 TV/m for larger numbers of electrons. 

Hence, radiation reaction is important for bunches containing > 10 6 electrons within a size <C A, that is, 
for bunches with a charge density close to solid state density. These electron bunches could only recently 



be extracted from atomic clusters by irradiation with an intense laser pulse. Of course, such dense elec- 
tron bunches quickly expand beyond R — A due to the strong Coulomb repulsion, so that they will scatter 
radiation coherently only for a brief period of time after creation. Therefore coherent Thomson scattering 
will not in itself be an efficient driving mechanism for electron bunch acceleration. However, it will be 
shown in this paper that even a short initial period of radiation reaction dominated interaction can have 
an observable effect. For larger electron bunches containing more charge, the radiation reaction effects 
are correspondingly stronger, and may be increased to such an extent that they can suppress the radial 
instability normally found in ponderomotive acceleration schemes, thereby enabling the full potential of 
laser- vacuum electron bunch acceleration to GcV energies. 



Collective radiation reaction of subwavelength electron bunches is, in some respects, the coherently 
enhanced version of the radiation reaction of a point charge. The latter has been a subject fundamental to 
the understanding of elementary par t icles and early electron theory and played a ro le in the development 



of quantum theory (jRohrlich . 



2007 



Lorcntz 



1916 



Abraham . 



1923 



Dirac . 



1938). It has always been 



a purely academic subject due to the smallness of the effect, h owever, and is only n ow coming within 



20041) . We suggest that 



experimental reach due to the availability of ultra-intense fields (jBulanov et al. 
an interesting alternative may now be offered by electron bunches of subwavelength size. As explained 
above, for bunches the enhanced radiation reaction effects are much stronger so that they are accessi- 
ble already using moderate laser intensities. Furthermore, also the characteristic time scale on which 
radiation reaction effects play a role (as measured by the quantity r e introduced below) is coherently pro- 
longed, bringing these effects into the experimental attosccond regime. Thus the coherent enhancement 
of radiation reaction by electron bunches yields interesting possibilities to investigate this fundamental 
subject by experimental techniques. 



In sections [2HH the bunch dynamics are analyzed treating the radiation reaction force in a perturba- 
tive way. In section [51 multiple-scale analysis and averaging over the laser optical time scale are applied 
to the equation of motion including radiation reaction, yielding a cycle-averaged description of the bunch 



dynamics. The results are applied to the case of a laser pulse in section [31 validating the averaged de- 
scription and our numerical calculations. Subsequently, in section 0] the theoretical results are placed in 
an experimental context by exploring the possibility of having subwavelength electron bunches to begin 
with, addressing both the time in which electron bunches Coulomb expand beyond wavelength and the 
production of dense bunches in laser-cluster interactions. Using accordingly realistic electron bunches, 
the observable effects of coherent radiation reaction on the bunch trajectory are calculated. Finally, in 
section[5]thc theoretical and numerical results are tentatively extrapolated to the non-perturbative regime 
of highly charged electron bunches. It is shown that ponderomotive acceleration of electron bunches by 
laser pulses may be stabilized by radiation reaction. 



2 Averaged radiation reaction force 



Suppose that an electron bunch of charge q = Ne, mass m — Nm e and radius Rq >C A oscillates 
due to the Lorentz force Fl of plane wave radiation propagating in the z-direction with electric field 
E = Eq sin (ut — kz) e x . The oscillatory quiver motion of the bunch (henceforth modeled by a point 
charge q) is a 'figure-of-eight '-cycle in the [x, z)-plane in the Lorentz frame in which the charge is at rest 



on average ([Landau &: Lifshitz 



20071) . In the course of each cycle, the accelerating charge continuously 
radiates and exchanges energy with its periodically changing Coulomb field, thereby experiencing an 
additional radiation reaction force Fr. The detailed descriptio n of Fr and th e associated equation of 



motion for a point charge have caused a long-standing debate dRohrlicb . 



equation of motion (|Lorentz 



1916 



Abraham, 



1923; 



Dirac 



1938) 



1997), in which the covariant 



(1) 



has played a central role. In eq. ([T]), dots denote differentiation with respect to proper time, is 
the electromagnetic field tensor of the incident radiation, r e = JV(2ro)/(3c) is the characteristic time for 
radiation reaction, and the metric signature (+,—,—,—) has been adopted. The first term in brackets 



equals the covariant Larmor power and represents the recoil of the charge caused by emission of radiation, 
that is, by the irreversible loss of four-momentum detaching from the charge and propagating towards in- 
finity. The second terms is referred to as the 'Schott term' and is associated with the reversible exchange 
of four-momentum with the Coulomb or velocity field, tha t is, with the part of the electromagne t ic fiel d 



1970; 



Her as fc O'Connell , 



2006). 



bound to the charge that does not give rise to radiation (jTeitelboim 
In the usual incoherent case N = 1, the bracketed terms can be neglected because of the smallness of r e . 
However, in the present case the radiation reaction parameter is TV times larger and / a m e iVm e , 
so that F£ oc N 2 . This shows again that the radiation reaction force is coherently enhanced and can 
become very relevant. Two regimes may be identified: the perturb ative regime in which lot e -C 1 and 



the radiation reaction dominated regime in which u)T e > 1 (see also (jBulanov et al 



20041 )). In the per- 



turbative regime, the Thomson cross section is by definition much smaller than the physical bunch size, 
and Fr <C Fl- This is the regime studied in this and the next two sections. In section [5J laser- vacuum 
acceleration in the radiation reaction dominated regime is considered. In passing, we note that eq. ([1]) is 
a fundamental equation of motion for a structureless point charge, so that one must be careful using it 
to describe a finite charge distribution, be it much smaller than the wavelength of the problem. In par- 
ticular, problems of interpretation occur when the mass equiva lent of the electr ostatic energy contained 



in the bunch becomes comparable to the rest mass m — Nm e fYaghjian 



20061 ). This is not the case as 



long as ujT e <C kRo, however, which is always satisfied in the perturbative regime. 



The two-dimensional 'figure-of-eight '-motion in the (%, z)-plane due to the linearly polarized incident 
radiation does not induce any y-component of acceleration in cq. |T]) , so that the motion of the charge 
remains two-dimensional throughout. Assuming Fr <C Fl, the radiation reaction force Fr affects the 
quiver motion negligibly on time scales ~ cu^ 1 . However, its cumulative effect over many cycles is to 
accelerate the charge, so that it gains a slowly varying average momentum p superimposed on the quiver 
motion. To study p, it is instructive to derive averaged equations of motion from eq. (fl]) by means 
of a multiple-scale expansion. For this purpose, define a fast dimensionless time scale T = uit and 
a slow dimensionless time scale eT = (wr e )T. The latter equals the time scale of radiative damping 



of a charged harmonic oscillator (Panof skv fc Phillips! . I1962D so that this is the time scale on which p 
will change appreciably. Hence, the dimensionless momentum (7, P x , P y , P z ) = P' 1 = p^/imc) may be 
written as the sum of a slowly varying part P M (eT) = (P M ) and a rapidly varying part P fi (eT,T) = 
p(°)^(eT, T) + eP < - 1 ' M (eT, T) + 0(e 2 ), where (•) denotes the average over an optical period. Incidentally, 
the quiver motion in a plane wave turns out to be more easily described in terms of the laser phase 



V 



T — kz rather than the time T explicitly (jLandau fc Lifshitz 



20071 ). Anticipating this fact by 



changing the fast coordinate from T to 77, time-differentiation takes the form 



y-±P»(eT,r,)=( n ° 







(2) 



where 7 2 = 1 + P% + Py + P| is the usual Lorentz factor. Substituting the sum Pf* = P P + in the 
space and time components of eq. ([T]), using ([2]), and collecting equal powers of e, yields the zeroth-order 
equations 



Pj°> 
p(o) = 7 (o) 



T 2 



—a cos 77; 

aP x cos 77 a 2 cos 2?7 
~ 1-Pz + 4(7-P 2 ) 
l + a 2 /2 + pt + pt, 



(3) 
(4) 

(5) 



where a = (ePg) / {m e cuj) is the dimensio nless field strength 



Eqs. (J3j{4j) give the well-known plane wave 



1996). Eq. ([5]) shows that the oscillating 



quiver motion with average momentum p ([Salamin fc Faisall 
charge b ehaves l ike a particle with effective mass my/1 + a 2 /2 due to the energy contained in the quiver 



motion ( Kibble 



19661 ). Expanding next the components of eq. ([TJ to first order in e, substituting the 



zeroth-order results Eqs. ([5][5]), and averaging over the fast time scale, one finally obtains 



dP x uJT e a r _ — ^ 2 Px 



dT 2 v ' ' 7 



tf-P*Y-£; (6) 



The averaged dimensionless position X = kx is calculated from 

d(X) - P 

Eqs. ([5][8]) now give the evolution of the averaged position X and momentum P of the charge in the 
plane wave. First, consider the special case P x = 0. Then Eqs. ([6][7]) reduce, after multiplication by 

m e cuj, to 

where 7 = €qcEq/2 is the intensity of the incident plane wave. This clearly shows that the averaged 
reaction force takes the classical form of a radiation pressure I/c on an effective cross section or, cor- 
rected by a Doppler factor j 2 (l — /3 Z ) 2 . This Doppler factor reflects the fact that plane wave appears 
Dopplcr-shiftcd in the frame of the charge, with a corresponding change in radiatio n pressure. Eq. (E l 
may also be obtained by Lorentz-transformation of the non-relativistic force (XtI/c (jLandau fc Lifshitz 



20071 ). In the more general case P x ^ 0, the direction of the radiation pressure is no longer along the 
propagation direction of the incident wave. This is because the secondary, scattered radiation produced 
by the charge is Lorentz boosted along _P, so that the recoil from this radiation has a net transverse 
component when P x ^ 0. The latter appears in Eqs. (|6][7]) as the additional factors involving P x on the 
very right. Effectively, in the transverse direction the averaged reaction force appears as a frictional force 
proportional to the velocity, which will be important bel ow. Effective forc es similar to Eqs. (J6]{TJ) have 



been obtained earlier for a one-dimensional plasma sheet (jirin et al. 



2001). 



A typical experiment in which the time-averaged forces exhibited by Eqs. ([5][7]) are important will 
involve the use of a high- intensity laser pulse of length t£ focused to a waist of size Wq. The laser 
pulse will appear to the electron bunch as a slowly modulated plane wave if ctl 3> A and Wo 3> fc _1 . 
The averaged reaction force will be as in Eqs. (j6][7]) with the dimensionless field strength a(X,T) now 
evolving according to the laser pulse shape. In this setting, however, the electron bunch will in addition 
experience a ponderomotive force due the the presence of large intensity gradients in the laser pulse. This 
is an additional average force which should be properly added to the right hand sides of Eqs. (|6][7|). It 
can be shown, by a multiple-scale analysis analogous to the one described above, that the ponderomotive 
force can be accurately described by 



(fc-iVHfflSpQ] 
47 



(10) 



if \/{cr L ) ~ (kwoy 1 < 1 (jQuesnel fc Mora 



19981 ). When adding Fp to the right hand sides of Eqs. 



([BUT]), in effect two separate multiple-scale analyses are combined. This is only correct provided that 
the small parameters used in both analyses are of the same order, which in this case means that also 
uJT e ~ X/(ct l ) — (kway 1 . 



Rather than using averaged equations, the detailed trajectory of the electron bunch may be calculated 
by direct integration of E q . ([T]). However, it is well-known that this equation yields unphysical runaway 



solutions (j Jackson . 



1999), which may be disposed of by substitutin g the applied Lorentz forc e p M ps F£ 



in the right hand side of Eq. ([I) . This yields the space component (jLandau fc Lifshitzl . 



2007) 



dp 
~dt 



dF L 7 3 fdv 

*y—jr - — I — x ( v x f l) 



dt 



dt 



This expression can also be derived from quantum mechanical arguments (lFo r d ii fc i O'Connell 



and has been put forward as the exact equation of motion rather than Eq. (fT) 



Rohrlich 



20011 ) 



(11) 



19911 ) 



3 Radiation reaction in a laser pulse 



In this section the results of the previous section will be applied to concrete examples of laser- vacuum 
experiments. First, however, it should be stressed that an accurate calculation of the detailed charge tra- 
jectory in a laser pulse requires the field of the pulse to be described in more detail than the usual paraxial 



appro ximations (|Hora et al. 



2000), even in the absence of radiation reaction forces. In ([Quesnel fc Mora . 
19981 ) , a laser field representation in Fresnel-type integrals is used to calculate detailed charge trajectories 
in a laser pulse, which are then compared with averaged trajectories calculated from the ponderomotive 
force Eq. (|10j) using the laser field in paraxial approximation. The agreement between the two results 
was found excellent except for highly relativistic in itial velocities. In this p aper we will adopt the same 



approach, using the detailed field representation of ([Quesnel fe; Mora 



19981 ) when evaluating the Lorentz 
force Fl in Eq. (|11[) . while taking the simpler paraxial approximation to evaluate the field intensity a 2 
in the averaged equations (|5l ll0[) . As a reference, both field representations are shown in the appendix. 



To validate the results of the previous section and to illustrate the effect of the effective radiation 
pressure produced by radiation reaction, we consider an intense laser pulse propagating in the z-direction 
and polarized in the ^-direction, which has wavelength A = 1 /im, peak intensity 7=1- 10 19 W/cm 2 and 
pulse length tl — 200 fs, focused to a relatively large waist of size wq = 10A. This pulse contains 2.3 J of 
energy and has a peak dimensionless field strength a max — 2.7. An electron bunch of radius J?o«) and 
radiation reaction parameter wr e = 0.03 is considered, corresponding to a bunch charge of q = 0.4 pC. 
For these parameters, A/(ctx) ~ (kwo)^ 1 ~ uir e so that the averaged description should be valid. The 
electron bunch is placed at rest at the initial position Xq = (xq,zq) = (A,0) prior to arrival of the laser 
pulse, deliberately at an off-axis position arbitrarily set to one wavelength, in order to study the radial 
acceleration as well. In this (unrealistic) first example, we disregard the Coulomb expansion of the bunch 
and assume that it retains its coherent Thomson cross section during the entire interaction, in order to 
study the basic equations governing the radiation reaction effects. 

Fig. [T] shows the trajectory of the bunch and its momentum as a function of time as it is overtaken by 



the laser pulse, calculated using both the averaged description Eqs. (fST fTTTj) and the detailed description Eq. 
(fTTj) . Clearly, the agreement between the averaged and the detailed description is excellent, validating 
both the multiple-scale analysis and our numerical calculation. The figure also shows the trajectory 
and the momentum in the absence of radiation reaction, that is, in case of conventional ponderomotive 
acceleration. In the latter case the charge is quickly expelled in the radial direction because of the large 
radial intensity gradient, which is a well-known instability t hat s o far has limited the applicability of 



ponderomotive acceleration schemes (jStupakov fc Zolotorev 



20011 ). The theoretical escape angle 9 for 



ponderomoti ve acceleration by a radia tion pulse of a charge initially at rest, is related to the final Lorentz 



factor 7/ by (jHartemann et al. 



19951) 



9 = arccos , 



Hf - 1 



(12) 



, 7/ + 1 

In the present example 7/ = 1.96, yielding a theoretical escape angle of 55.3° in excellent agreement 
with the numerically calculated result. In contrast to the purely ponderomotive case, in the case includ- 
ing coherent Thomson scattering the acceleration in the radial direction is markedly suppressed. This 
stabilizing effect is due to the average frictional force shown by Eq. ([6]), which opposes the tendency to 
radial acceleration. Meanwhile, the charge stays in the beam longer and is accelerated in the forward 
direction for a longer time by both the radiation reaction pressure and ponderomotive forces, leading to 
a higher final energy and a lower escape angle. 



4 Realistic electron bunches 

We now consider the more realistic case of an electron bunch that expands due to Coulomb repulsion. 
The time after which the Thomson cross section fails to be coherent, that is, the time it takes the bunch 
to expand to R = A, depends on the charge and the initial size of the bunch. To obtain realistic values for 
these initial parameters, we first make some estimates concerning the method to produce dense electron 
bunches mentioned in the introduction. 



Rather than irradiating a pre-existing electron bunch, an electron bunch may be created in practice by 
irradiating a sub- wavelength atomic cluster. When the laser pulse is sufficiently strong, the electron bunch 
will be emitted from the cluster somewhere in the leading edge of the pulse, after which it is available to 
be accelerated by the remainder of the pulse as has be en analyzed in the prev i ous section. Many d i fferent 



mechanisms play a role in laser-cluster interaction (jSaalmann et 



2006 



Krainov &: Smirnov 



2002) 



and it is beyond the scope of this paper to a nalyze this interact ion in detail. Instead, we use a strongly 



simplified model of the cluster suggested by (jParks et al. 



20011 ) to obtain order-of-magnitude estimates 



for the initial parameters of the electron bunch. Accordingly, at some point in the leading edge of the pulse 
the field strength becomes such that all atoms in the cluster are ionized almost instantaneously, changing 
the neutral cluster into a dense plasma ball. Subsequently, this plasma ball is modeled as a spherical rigid 
electron bunch, interpenetrating with and moving through a practically immobile oppositely charged ion 
bunch, under the action of both the driving laser electric field and the restoring Coulomb force of the 
ion cloud. As the rising edge of the laser pulse advances, the electron bunch oscillates around the ion 
bunch center with increasing amplitude. At a critical laser field strength the electron bunch breaks free 
and escapes the ion bunch within a fraction of an optical period, yielding the free, dense electron bunch 
we will now use as an input for our calculation. 

The critical dimen sionless field strength at which the electron bunch is liberated may be estimated as 



( Parks et al 



2001) 



&crit 



(13) 



where u> p = y/'SNe 2 / (47reom e i?Q) is the plasma frequency of the ionized cluster consisting of N atoms. 
The electron bunch le aves the ion bunch at the maximum of an optical cycle. The time this takes is of 
the order t csc ~ a;" 1 ( Parks et al. 



2001), so that the escaping electron bunch will have a velocity of 



approximately 



(14) 



in the direction of laser polarization. Incidentally, note that the coherent radiation reaction parameter 
can be expressed in terms of the cluster parameters as ur e — (2/9)(tJ p /w) 2 (fci? ) 3 , by which the validity 
condition wr e <C kR of section [5] automatically implies that /3 esc < 1. After leaving the cluster, the 
electron bunch will start to expand. It can be shown straightf orwardly that, af ter a short period of slow 



expansion, the expansion rate increases to the constant value (jAmmosov 



19911) 



dR 
~dt 



1 + P 2 



(15) 



In Eq. (I15p . it has been used that the expanding electron bunch will be in quiver motion immediately 
after leaving the cluster, so that it will move with an average Lorentz factor according to Eq. (|S|). This 



reduces the expansion rate by a factor 7 = \Jl + P 2 SC + a 2 Yit /2. 



As an example we take an i?o = 16 nm cluster with an electron density of 5-10 m , which can be rou- 



tinely made using a supersonic gas jet (jKrainov fc Smirnov 



2002) , irradiated by the same laser pulse as in 



the previous example. With these numbers the other parameters are fci?o = 0.1, u> p /u> = 6.6, wr e = 0.01, 
and a crit = 0.68 and (3 CSC — 0.66 according to Eqs. (|13H14j) The last two numbers roughly compare t o 



2010: 



Fukuda et al 



20071) 



the field strengths used and electron velocities obtained in (jLisevkina et al 
Assuming the expansion rate Eq. fjl 5f) . the bunch grows larger than the wavelength after a time 8.3w _1 , 
that is, after no more than 1.3 optical cycles following the moment of escape from the cluster. Thus 
coherent Thomson scattering will only occur for a brief initial period. Nevertheless, even such a brief 
period can have an observable effect on the trajectory of the electron bunch in interaction with the laser 
pulse. 

Fig. [2] shows the trajectory of the bunch in interaction with the laser pulse, both for the case without 
the radiation reaction force (wr e = 0) and for the case including the radiation reaction force (wr e = 0.01) 
during the initial period < cut < 8.3. In the latter case, the decrease in coherent radiation reaction 
was modeled by a sudden switch-off, setting wr e = at the somewhat arbitrary point cut = 8.3. Again, 
the results of Fig. [2] were calculated both using the averaged description Eqs. (j5lTlO|) and the detailed 
description Eq. (JTTJ) . In both cases the initial position was Xq = (— A, 0) and the initial velocity was 



(3 esc ce x . The initial position of the laser pulse along the z-axis was chosen such that a = a cr i t at xq, and 
the optical phase offset </> was chosen such that the optical cycle was at its maximum at t = 0. This 
time, the escape angle for the case without radiation reaction does no longer agree with Eq. (fT2|) or its 



generalization for a nonzero initial velocity (jHartemann et 



19951 ). This may be caused by the sudden 



injection into the radiation field at t = 0, which is not considered in the derivation of Eq. (fT2j) . The good 
agreement between the averaged and detailed descriptions is remarkable here, considering the fact that 
the concept of averaging the equation of motion over the optical time scale is not a priori valid when 
applied during the initial period, which is as short as the optical period itself. 

From the figure it is clear that the initial period of radiation reaction still reduces the escape angle 
observably, be it much less than in the previous example. To put the 1.3 degrees of reduction in angle 
in perspective, the figure also shows the error margins around the trajectory without radiation reaction 
which would be caused by an uncertainty of magnitude A in the initial position Xq. The deviation caused 
by radiation reaction is clearly outside these error margins. 



As another example we show how the initial period of coherent Thomson scattering may even change 
the behavior of the bunch completely. Fig. [3] shows the trajectory of the bunch, using the same laser 
parameters and initial conditions as in the previous example, only now with the initial position changed 
to Xq = (— 7.1A, 0). This time the effect of radiation reaction is such that the bunch is deflected by the 
laser pulse in the negative x-direction instead of the positive x-direction. Under the initial conditions of 
this example, the bunch is produced much further from the laser beam axis and is therefore decelerated 
in the negative x-direction by the radial ponderomotive force of Eq. (|TU| for a longer time. Without 
the action of the radiation reaction force, the bunch just makes it across the laser axis (where the 
ponderomotive potential is maximum) , after which it is accelerated in the positive x-direction. Including 
radiation reaction, however, in the initial stage the bunch is additionally decelerated by the frictional 
force of Eq. ([5]). Furthermore it is additionally accelerated in the positive z-direction by the radiation 
pressure of Eq. (J7J , so that it keeps up with the laser pulse and is under the action of the decelerating 
radial ponderomotive force for a somewhat longer time. Both effects are only small, but in this case just 



enough to prevent the bunch from passing the beam axis. It is to be noted that in this example whether 
or not the bunch will pass the beam axis is very sensitive to the initial parameters. This is an additional 
reason why the averaged description cannot be used here to reproduce the trajectories of Fig. ([3]), since 
the small differences between the two descriptions in the initial part of the trajectory affect the final 
behavior strongly. Still, it is clear that radiation reaction can play an important role in laser-vacuum 
experiments. 

5 Pre-acceleration using radiation reaction 

In the previous section, relatively small electron bunches were considered leading to modest values for 
the bunch charge and for the radiation reaction parameter (wr e <C 1), thereby staying well inside the 
perturbative regime studied in sections [2H1 Much stronger radiation reaction effects could be expected 
when scaling the system to larger bunches. In the introduction, the magnitude of the effective radiation 
reaction force was estimated assuming a wavelength of I /im. Taking now, for instance, a CO2 laser with 
a wavelength of 10 (im, ten times larger bunches still Thomson scatter coherently. Since the coherent 
Thomson cross section scales as <tt oc TV 2 oc i? 6 , in that case an effective force of Ft = (JtI/c ~ kN 
rather than mN may be expected, equivalent to accelerating fields approaching Fp j (Ne) ~ PV/m rather 
than TV/m. Thus, the radiation reaction force in this case completely dominates the bunch dynamics. 
But as before, the action of the radiation reaction force is restricted to the initial period following the 
creation of the electron bunch, while the bunch is still smaller than the wavelength. Therefore the large 
accelerating fields just mentioned are not effective as a driving acceleration mechanism. However, the 
initial radiation reaction dominated phase pre-accelerates and redirects the velocity of the bunch, which 
may serve to stabilize subsequent ponderomotive acceleration. 

As an example of such acceleration using radiation reaction, consider the same laser pulse as used in the 
previous examples, but scaled ten times in all directions so that A = 10 /im, tx = 2 ps and wo = 10A. 
This CO2 laser pulse contains 2.3 kJ of energy and has a peak dimensionless field strength a max = 27, 
which is presently available. The increased wavelength allows an electron bunch of initial size Rq = 160 
nm containing N = 8.5 • 10 8 electrons, resulting as before in kRo = 0.1, but this time in a much bigger 



value 0JT e = 1. The non-relativistic cluster model of ( Parks et al 



2001) does not apply for such a large 



amount of ch arge. We assume a rel ativistic initial electron bunch velocity of /3q = 0.9 instead, as is also 



indicated in ( Lisevkina et al. 



2010h . Substituting P csc = (3 / yj 1 — /3g and a cr i t = a max in Eq. (|T5l) . the 



resulting expansion rate is such that the bunch stays smaller than the wavelength for about one optical 
cycle, just as in the previous examples. 

Fig. [4] shows the trajectory and momentum of this electron bunch in the laser pulse just described. The 
initial bunch position was set to a;o = (A,0), the initial position of the laser pulse along the z-axis was 
chosen centered around x-axis so that the bunch started its movement at the pulse maximum, and the 
optical phase offset 0o was chosen such that the optical cycle was at its maximum at t = and Xq. From 
the figure the effect of the initial period of radiation reaction is evident: the large radial velocity of the 
bunch is suppressed immediately preventing an early escape from the beam waist, after which the bunch 
is strongly accelerated in the direction of propagation of the laser pulse. The inset of the p z -panel of 
Fig. 4 shows the long-term evolution of the momentum in this direction. The bunch is first accelerated 
by the leading edge of the laser pulse and then decelerated by the trailing edge, as is characteristic for 
laser- vacuum acceleration. The energy of the bunch is increased to 7 ~ 600, corresponding to an energy 
of 0.3 GeV per electron. 

In the case without radiation reaction, the bunch is quickly expelled from the laser beam. Consequently, 
the energy gained is much l ess. The cusp-like features in the corresp onding p z - and p x -plots of Fig. U 



agree with previous results (jDodin fe Fisch 
motive scattering. 



2003 



Galkin et al 



20091) in comparable cases of pondero- 



Returning to the case including radiation reaction, it may seem strange that the bunch starts, at 
t = 0, with a downward radial acceleration while the laser electric field points upward at that moment. 
This, however, is a consequence of the large value of ujt e as can be inferred from Eq. (1111) . Namely, at 
t = the optical field at the bunch position is maximum, so that 



qE 



e x m e cuja ma , x ; 



qcB 
dF L 

dt 



G y 771 e CCU CI m ax , 

dv 

q Tt xB ' 



while p z = 0. Substituting the above in the x- and z-components of Eq. (jlip . writing dv / dt in terms 
of the momentum p = ~/mv, and rearranging, one finds in dimensionless form the x-component 



dP x 
dT 



1 - WTeflmax 7qA) 



t=0 1 + 7o( WT ea m a X )^ 

where 70 = (1 — /3q) -1 ^ 2 - In the example of Figure 21 wr e a max/ 9 ^ 1) so that Eq. (flU]) reduces to 



(16) 



(17) 



In the inset of the p^-panel of Fig. IU this initial acceleration is indicated by the sloped line. Eq. (p~6|) 
suggests that the initial velocity of the bunch can be redirected by means of a short period of radiation 
reaction choosing favorable parameters a max and wr e , and may be used to efficiently accelerate the bunch 
into the direction of propagation of the laser such as has been done in the example of Fig. |U 
It should be noted that the very large radiation reaction parameter used in the last example poses some 
challenges to the theory of this paper. In particular, here wr e ^> kRo, which makes the applicability of 
the point charge equation of motion (fT]) to an extended electron bunch debatable. Furthermore, Eq. ([Pi]) 
has been proposed on the o ne hand as an approximat ion to the classical equation of motion (TTJ) based 



on the assumption ujt e <C 1 (jLandau fc Lifshitz 



2007) , and on the other hand as the exact equation of 



20011 ) . Interestingly, the experimental availability 



motion replacing Eq. ((T|) for all values of uiT e (jRohrlichl 
of charged objects for which wr e ~ 1 may now give some insights about which approach is right, and may 
offer further opportunities to experimentally test radiation reaction theories. 



6 Conclusion 



The current developments in laser technology make it possible to obtain subwavelength electron bunches of 
very high charge density. With that, a qualitatively new regime is accessed in which coherently enhanced 
radiation reaction effects become significant. In this paper, we have analyzed some of these effects in the 
context of laser- vacuum experiments. It has been shown that the radiation reaction force affects the bunch 
dynamics notably, even if the radiation reaction can still be treated as a small perturbation on the optical 
time scale. Considering larger bunches containing more charge, we demonstrated that the radiation 
reaction effects may even become strong enough to be exploited in effective bunch acceleration schemes, 
although this is also where the theory needs to be further developed. It is clear that the coherently 
enhanced radiation reaction of high-density electron bunches offers interesting new possibilities, both as 
a technological tool in the development of novel acceleration schemes, and as an experimentally accessible 
system to study the fundamental topic of radiation reaction. 

Appendix A: Laser field representations 



In th is paper, when evaluating the Lorentz force Fl in Eq. (JTJ), the following fields are used ([Quesnel fc Mora 
199a): 



cBy = ^^U-T^ + Va); (20) 



= I^^J 4 ; (19) 
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(j> — uit ~ kz\/l — b 2 + 4>o, 

where 4>o is the phase offset of the optical field. These expressions describe a Gaussian beam which 
propagates along the z-direction, is polarized in the ^-direction, and is focused to a waist of size wq. To 
obtain a laser pulse of length tl instead of a full beam, the above fields are multiplied by an envelope 
function cos 2 (wrj / (2ujtl)) ■ To evaluate the field intensity a 2 in the averaged equations (l5HT0|). the simpler 
paraxial approximation is used according to 



^^ e -^t(M 2 cos 2 ( El V 
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Figure 1: a) Trajectory of an electron bunch accelerated by a laser pulse with parameters A = 1 /Ltm, 
1=1- 10 19 W/cm 2 , tl = 200 fs, and w = 10A; b) Longitudinal momentum of the bunch as a function of 
time; and c) Transverse momentum of the bunch as a function of time, calculated for the case including 
radiation reaction (Thomson + ponderomotive scattering, ujT e = 0.03) and for the case without radiation 
reaction (ponderomotive scattering, lot b = 0). The solid lines have been calculated using Eq. (fTTj) with 
Fi = q(E + v x B), taking for E and B Eqs. (|18H21[) . The dashed lines have been calculated using Eqs. 
([5l fTU|) . taking a 2 according to Eq. ([22]) . The initial position was Xq = (A,0) and the initial velocity was 
zero. 



Figure 2: Trajectory of an electron bunch accelerated by a laser pulse, calculated for the case including 
radiation reaction (Thomson + ponderomotive scattering, wr e = 0.01) during the initial period < ut < 
8.3, and for the case without radiation reaction (ponderomotive scattering, lot c — 0). See Fig. (fTJ) for 
details about the calculation and laser parameters. The dash-dotted lines show the error margins for 
the case wr e = corresponding to an hypothetical uncertainty of magnitude A in the initial position, 
calculated using the averaged description. The initial position was x = (—A, 0) and the initial velocity 
was according to Eq. (fT4|) . 



200 




200' 

50 100 150 200 

kz 



Figure 3: See Fig. ((2J, with the initial position changed to x = (— 7.1A, 0). 



Figure 4: a) Trajectory of an electron bunch accelerated by a laser pulse with parameters A — 10 ^tm, 
1=1- 10 19 W/cm 2 , T£ = 2 ps, and u>o = 10A; b) Longitudinal momentum of the bunch as a function of 
time; and c) Transverse momentum of the bunch as a function of time, calculated for the case including 
radiation reaction (Thomson + ponderomotive scattering, wr e = 1) during the initial period < cut < 2tt 
and for the case without radiation reaction (ponderomotive scattering, u)T e = 0). The solid lines have 
been calculated using Eq. (JTTJ) with F L = q(E + v X B), taking for E and B Eqs. (|18H21j) . The inset 
of the middle panel shows the long-term behavior of the longitudinal momentum for the case including 
radiation reaction. The inset of the lower panel is a zoom-in on the first optical period; the slope of the 
dashed line is according to Eq. (flTj) . The initial position was Xq = (A,0) and the initial velocity was 
Po = 0.9e x . 



